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Introduction 

The motion of a charged particle in electromagnetic fields 
of various geometric configurations and arising from a variety of 
sources is of intrinsic interest in electromagnetic theory. The 
results of the analysis of the motions generated by the field- 
charge interaction are in many cases also pertinent to the 
discussion of applied problems in plasma physics and astrophysics. 
The particular configuration consisting of a plane wave propagating 
in the presence of a : static uniform magnetic field whose direction 
is parallel to the wave normal is examined in this report. The 
analysis presented here is developed along lines different from 
those pursued in a previous study of the same configuration by 
Roberts and Buchshaum [1964]. It should be remarked that the 
problem under consideration is treated within the context of 
classical electromagnetic theory. The corresponding quantum 
mechanical problem - but neglecting the quantization of the field - 
has been discussed by Redmond [1965] who solved the appropriate 
Klein-Gordon and Dirac equations. In neither of the references 
cited above Is the radiative reaction force taken into account. 
While from a physical point of view a knowledge of the influence 



of the radiative reaction on the motion of the charge would be 
very useful, the analytical difficulties associated with the 
Lorentz-Dirac equation are frequently sufficiently formidable to 
preclude the obtaining of the desired information. In fact it 
was hoped originally, on the basis of some success with a simpler 
field configuration, Mitchell et al. [1971], to include this 
radiative reaction in the following analysis, but to date mathe- 
matical difficulties have proved insurmountable in the search 
for a purely analytical solution. A numerical solution - at least 
to the approximate Lorentz-Dirac equation - could, of course, be 
obtained in a relatively straightforward way. 

As a final comment prior to developing the analysis, it may 
be pointed out that an approach to this problem through the 
Hamilton-Jacobi formulation, although not adopted here, is an 
alternative possibility. Such an approach would naturally have 
much in common with the quantum mechanical formulation. The basic 
non-linearity of the Hamilton-Jacobi equation, especially in its 
relativistic form introduces major difficulties, however; see 
e.g., Corben and Stehle [1965], Landau and Lifshitz [1959]. 
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Motion in Moving Frame of Reference 


The direction of propagation of the wave, and consequently 
the uniform magnetic field direction, is taken to be along the 
positive z-axis of the usual cartesian xyz triad. The field can 
then, without significant restriction, be represented as the 
superposition of a circularly polarized monochromatic wave and a 
uniform field, H , in the form 


E = E sin (ut - az) 

X 


E = -E cos (tot - az) 


E = 0 
z 


H = H cos (tot - az) 

X 


H = H sin (tot - az) 

y 


H = H 
z o 


(1) 


all quantities being measured relative to a reference frame in 
which the electrical conductivity is zero and the constitutive and 
field equations are the usual 


... ( 2 ) 


eE ; 

V-H = 0 ; 

VxE + vi 

= o 
3t 

UH ; 

o 

II 

• 

> 

VxH - e 

= 0 
3 1 


expressed in rationalized m.k.s. units. It follows ,. Stratton [1941], 
that the phase velocity, v^, of the wave .is given by 


v = to/a = 1/v/yF = E/uH 


(3) 


and that the index of refraction, n, of the medium in which the 
motion takes place satisfies the relations 

n = c/v^ = c/to = cB/E (4) 

where the speed of light in vacuo is represented by c. Subsequently, 


4 


unless explicitly stated otherwise, it is assumed that n > 1. The 
exception to this inequality will be the case in which n = 1, i.e., 
the case in which the motion takes place in a vacuum. 

The equation of motion of the charge, q, is 


—2. = q[E + vxB + vxB ] (5) 

at ^ o 

in which the momentum of the charge is 

->• / 7 2 2 

p = mv//l — v /c z (6) 

m, being the rest mass. 

The analysis of the motion can be simplified considerably by 
carrying out a Loren tz Transformation to a frame moving in the 
positive z-direction with velocity v^. Relative to this frame the 
particle is subjected to the influence of a magnetic field alone, 
the electric field of the wave being entirely eliminated. In 
fact, after some reduction one finds the magnetic field in the 
moving frame to .be 


i, yH 

B' =1 •=— COS 


V 


co . t uH . u ' i ■ u 
z ' - j ■=— si-n - ■ z 1 + kpH^ 


1*1 


v i*i 


and the equation of motion reduced to 


dp! _ x g, 

dtT - q v x B 


(7) 


( 8 ) 


In equations (7) and (8) the superscript prime has been introduced 

to denote quantities measured relative to the moving reference 

frame, the unit vectors i , j, & are along the x' , y', z' axes, 

2 2 ” 1/2 

respectively, and the constant factor, represents (1 - v^/c ) 

Equation (8) in cartesian component form is 
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dP ' J , 

dtT = °2 dt 7 " + a l d tT Sin 0 3 2 ’ 


dz 


^ p y dz 1 „ , dx' 

dt^ =a ldF cos a 3 Z °2 dF 


and 


dp' 

dt' 


dx • , dy* 

J 1 dt^ Sin a 3 2 “ °1 dt 7 * 


dZ_ cos a z • 
1 dt' 3 


the three constants 


a 1 = qpH/yj^ , a 2 = qyH Q . , . c 3 = w/y 1 v 1 


(9) 

( 10 ) 


(ID 


( 12 ) 


being introduced for conciseness. 

The first two component equations can be integrated immediately 
to give 

p x = a 2^' ” a 4 cos a 3 2 ' + C 1 (13) 

and 

Py = - a 2 x ’ + a 4 sin a 3 z ' + C 2 (14) 

respectively, where and C 2 are constants of integration and the 
constant = o^/a^. 

The energy of the charge is 


E' = me 2 /A ~ . (^-) 2 ' 

and so 

ML - v . ; ML 

dt' v dt' 


(15) 


(16) 


From equations (9), (10), (11), and (16) it follows that 


dS 1 

dt 7 - 


0 


(17) 
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and therefore the energy E' is a constant. This, of course, is to 
be expected because the motion is taking place in a pure magnetic 
field. Accordingly, equations (13) and (14) can be written 

my' = 0 2 y '“ oj cos a 3 z ' + C 1 (18) 

and 

my ' =-0 2 x'+— sin0 3 z'+C 2 (19) 

_1 

v i 2 "" 2 

in which the constant y ' = (1 - (-^-) ) .' Although the third compo- 

nent of the equation of motion cannot be integrated as directly as 
the first two, the explicit spatial deperidence of p' can be found 
as follows. On eliminating the terms dependent on z' from the 
right hand side of equation (11) one finds after some development 
an integrable equation which leads to ... 

p z = " I a 2° 3 (x ' 2 + y ' 2) + a 3 (C 2 X ' " C l y,) + C 3 (20) 

and consequently 

my ' ~r = - j o 2 a 3 (x’ 2 + y' 2 ) + ^(C^x- - C^') + C 3 (21) 

in which is a constant of integration. ■ 

Alternatively, the temporal and spatial and dependence of 
p' can be determined from the general energy -momentum relationship 

E' = (p'c) 2 + m 2 c^ (22) 

The trajectory of the charge in the moving frame of reference is 
given in parametric form by equations (18), (19), and (21), or in 
simultaneous form by 
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dx 1 

0 2 y* - a 4 cos a 3 z 1 + 


dy | 

-c^x 1 + sin a^z 1 + C 2 


dz 1 

-a 2 a 3 (x' 2 + y ,2 )/2 + a 3 (C 2 x' - C^y’) + C 3 . 


(23) 


Without attempting an explicit determination of the trajectory it 
can be asserted on the basis of equations (13) , (14) , (20) , and 

(15) that the motion takes place on the surface 

/ 

~ 0 2° 3 2 2 1 2 

— 2 — ( x ' + -y ■ ) - a 3 ( C 2 X ' “ c iy’) “ C 3 J 

2 2 

■+ (c^y 1 * ~ °4 cos a 3 Z ' + C l^ + ^ a 2 X ' ~ a 4 S ^ n a 3 Z ' “ C 2 > 

= constant (24) 

where the constant is determined by the initial conditions. 



8 


Motion in Fixed Frame of Reference 

The application of the Lorentz Transformation to the momen- 
tum four-vector 


M ' = {my 


. i 51> 

c ' 1 2 s 

c 


(25) 


produces the equations which describe the motion in the fixed frame, 
After some algebraic manipulation one finds the momentum components 
to be 

P x = a 2 Y - cr^ cos (ait - az) + (26) 


p = ~<^ 2 X ~ a 4 s i n (wt- - az) + C 2 


(27) 


and . 


p z - -V 1 o 2 o 3 (x 2 + Y 2 )/2 + r 1 a 3 (C 2 x - c x y), + YjCj + y^ ^ 

(28) 

The energy of the motion relative to the fixed frame can be written 
in the form 

. E = y ± E' + y iVi p; 


= y ± E' - 'Y 1 v 1 a 2 a 3 (x2 + Y 2 ) /2 + Y x ~ CqY) + Y i v l C 3 

(29) 

It is worthy of note that equations (28) and (29) can be combined 
to give 


E' = Y 1 (E - v x p z ) 


(30) 


Thus the quantity E - v^P z remains constant throughout the motion. 
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An additional constant of the motion can be obtained by combining 
equations (26) and (27) in the form 


(?,. - a 2 y - C 1 )“ + (p y + c 2 x - C 2 ) = c 4 


(31) 


A more direct link between the relationship expressed in 
equation (30) and the equation of motion (5) may be established as 
follows. The rate of change of the energy E is 


C-2 ±. CLjD 


at 


- v 


dt 


(32) 


wnxcn 


= qv • 2 


(33) 


on using the equation of motion (5).. However, for a plane wave 


2 = Vl B x k 


and again from equation (5) 


■ j-* dt 

ClP “2 r>- , ^ » 

K * dt = “dt = * (v * B) 


(34) 


(35) 


Therefore , 




dp 


= v 


1 dt 


(36) 


ana 


-2 - V- (p - p ) 

O X £ 2 ‘ Z 

O 


(37) 


in which the subscript zero- denotes values at some particular time 
t o - It is now seen that the constant 2 ' introduced in equation (15) 
satisfies the equation 


S' = YJ_ (Z - v x p ) 

o 


(3S) 
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In conclusion one notices that a simultaneous set of 
differential equations which determine the trajectory relative to 
the fixed frame can be written by combining equations (26), (27), 

and (28) . In contradistinction to equation (23) , however, the 
formulation under discussion would explicitly involve the time. 
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Motion in a Vacuum: The Case n = 1 

The analysis contained in the preceding pages has been 
developed subject to the restriction that the index of refraction, 
n, of the medium in which the motion takes place is greater than 
unity. However, the case in which n = 1 is also of considerable 
interest. It cannot be treated by the Lorentz Transformation 
already used because the phase velocity of the wave is v^ = c 
when n = 1, and hence y^ -*■ °° rendering the transformation singular. 
One observes, nevertheless that the expressions (26), (27), and 

(28) for the components of the linear momentum are independent* 
of the factor y It can be shown, in fact, by direct substitution 
in the equation of motion (5) that the three relations (26), (27), 
and (28) are valid for any value of n. Accordingly, the constants 
of the motion already established exist independently of any 
restriction on the value of the index of refraction. For the 
case n = 1, the energy relationship (37) takes the form 

E = E + c(p z - p z ) (39) 

o 

On using the general energy-momentum relationship, one finds that 

p z - \ <p z -V c) ' 1{(p z - V c)2 ■ m2c2 

o o 

2 

- [a 2 Y ~ cos (wt - az) + C^] 

2 

- [a 2 x + a 4 sin (wt - az) - C 2 ] } (40) 

* 

Equation (28) exhibits a superficial dependence on y^- 
removed upon recourse to equation (12) . 


It can be 
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Equation (40) is equivalent to equation (28) which, although derived 
for n .> 1, is now seen to be valid for arbitrary n. 

In addition equation ( 40) when used in conjunction with 
equation (39) provides the explicit dependence of the energy E on 
the position of the charge and on time. This functional dependence, 
as in the case n > 1, can be expressed in the form 2?(x, y, cot - az) 
without difficulty. 
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Motion in Momentum Space 


Considerable insight into the particle’s motion can be 
obtained by examining the dynamics in a "momentum- space " in which 
the momentum components are taken to be the rectangular cartesian 
coordinates of a point. In particular, the difference in the 
behavior of the charge in the two cases n > 1 and n = 1 is quickly 
made apparent by the formulation in momentum-space. For this 
reason, and also for mathematical convenience, the two cases referred 
to are once again treated separately. 


1 . Index of Re f Taction n >_ l_ 

On squaring equation (37), which actually expresses the 
constancy throughout the motion of a certain function of the momen- 
tum components, it is seen that the charge is confined in momentum- 
space to the surface defined by 

v n 


2 2 1 2 v 1 22 
p + p + —rr p - 2 —75- ( E - V.p )p + in c 

*x y 2 c 2 o ± Z Q z 


4 (£ o - v l p z > 2 ■ 0 
c o 


(41) 


This equation represents a prolate ellipsoid whose axis of symmetry 

coincides with p -axis. The semi -major axis, a, of this spheroid 

z 


IS 


a = (P 2 - P z , )/2 - Yl ((S 0 - v lP ) 2 (-^) 2 - m 2 o 2 J 

max mm o 


1/2 


where the maximum and minimum values of p are, respectively 

z 


max 


v- 


■£+ [<|-r - m-c*] 


2 2 


1/2 


(43) 
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ram 


n 


{ £1 . - [(£1)2 - m 2 c 2 i 1/2 } 

L c c c 


( 44 ) 


Owing to the closed nature of the surface of the ellipsoid 
it is immediately clear that the energy of the charge remains 
finite -and bounded at all times- In fact, it is possible, although 
the results will not be reproduced here, to calculate, explicit 
upper and lower bounds on the energy. 


2. Motion in a_ Vacuum ; Index of Refraction n_ = 1_. . 

In this case the appropriate point to start the discussion 
form is equation (39) which leads, on squaring it, to 


P + • p 

-x 


E 

2( ~r " 


)p. 


, 2 2 
+ m c 


E 

- <-f - 


> 2 

P Z } = 

o 


(45) 


The significant difference between equations (41) and (45) is the 
■ absence of a quadratic term in p z in the latter. This 
results in an opening up of the ellipsoid into a paraboloid, its 
surface being a surface of revolution about the p z ~axis. 
Consequently, it is no longer possible to establish an upper 
bound on the particle's energy which may therefore increase with- 
out limit as the motion develops. Some aspects of these results 
are discussed in a different context by Roberts and Buchsbaum 
[1964] . 

As a final comment, it is remarked that if the index of 

refraction n is less than unity, the corresponding surface in 

momentum- space is determined by an equation formally similar to 

equation (41). However, the condition n < 1 implies, that the 
2 

coefficient y-j_ < 0 and hence the surface is a. hyperboloid of 



r <3v< 


°l Ut 


l o n. 


u Ppe r 


As 


^°u nd 


in t hi 


on 


the 




case 

9y. 


n = 
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